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Numerically  Efficient  Use  of  Frequency  Domain 
MoM  Codes  For  Wideband  Radiation  and 
Scattering  Problems  -  An  Application: 
Scattering  From  a  Microstrip  Patch 


1.  miRODUcnoN 

with  the  present  Interest  in  electromagnetic  pulses  the  question  arises  whether  existing 
frequency-domain  computer  codes  can  be  used  to  provide  the  temporal  response  of  antennas  and 
scatterers. 

These  codes,  which  usually  are  based  on  the  Method  of  Moments  (MoM),  are  highly  accurate 
and  also  represent  a  large  investment  in  both  flnancial  resources  and  development  time.  Of 
course,  since  time  and  frequency  domain  solutions  are  related  by  the  Fourier  transform,  these 
codes  in  principle  could  be  used  to  numerically  evaluate  the  spectrum  of  the  impulse  response 
and.  after  multiplication  with  the  signal  spectrum,  an  inverse  Fourier  transform  would  lead  to  the 
desired  time  domain  response.  However,  the  impedance  or  scattering  behavior  of  many  objects  is 
not  a  smooth  function  of  frequency,  especially  at  frequencies  near  resonances  of  the  structure, 
and  to  capture  these  rapid  variations  requires  a  very  dense  numerical  sampling  of  the  spectrum. 
This  represents  a  large,  and  sometimes  prohibitively  large,  computational  effort. 

Thus,  there  is  a  need  for  a  more  efficient  approach  than  a  straightforward  numerical  Fourier 
transform. 

The  approach  proposed  here  is  based  on  the  observation  that  the  late  time  impulse  response 
of  many  3-dimenslonal  structures  is  often  well  represented  by  a  set  of  damped  sinusoids  |1]  that 
appear  as  singularities  (poles)  in  the  complex  frequency  plane.  This  suggests  modeling  the 
response  by  a  small  set  of  such  complex  natural  resonances  and  associated  natural  modes,  rather 
than  by  a  purely  numerical  spectrum  evaluation.  The  approach  can  be  viewed  as  "smart" 
interpolation  since  the  physically  motivated  Interpolation  functions  contribute  additional 
information.  The  method  has  also  been  called  model-based  parameter  estimation  [2],  a  term 
borrowed  from  signal  estimation  theory. 
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The  erall  approach  involves  the  following  three  general  steps.  First,  an  appropriate  number 
of  complex  natural  resonances  for  the  frequency  range  of  Interest  and  the  corresponding  natural 
modes  are  determined.  This  is  done  with  the  aid  of  a  frequency  domain  moment  method  (MoM) 
code  assumed  available  for  the  structure  under  consideration.  Second,  the  induced  currents  on 
the  structure  are  expanded  in  terms  of  the  natural  resonances  and  modes  with  unknown 
amplitudes.  This  is  the  so-called  singularity  expansion  method  (SEM)  current  estimate.  The 
determination  of  the  mode  amplitudes  is  typically  the  most  difficult  step  of  the  overall  process  (IJ, 
[3].  {This  is  where  people  usually  give  up',  to  quote  Dr  D  Girl,  an  old  master  of  the  art  of  the 
SEM).  Presently,  a  new  approach  is  Introduced  to  determine  the  amplitudes  that  uses  least- 
squares  matching  of  the  current  estimate  to  known  or  'trusted  reference'  currents  at  several 
frequencies.  The  reference  currents  are  produced  with  the  existing  MoM  code.  Finally,  the 
electromagnetic  characteristics  of  interest  such  as  the  input  Impedance  or  radar  cross  section 
(RCS)  are  calculated  from  the  SEM  current  estimate.  This  can  be  done  rapidly  over  a  wide 
frequency  range,  due  to  the  simple  form  of  the  current  estimate. 

One  of  the  advantages  of  our  method  includes  the  possibility  of  using  existing  frequency 
domain  MoM  codes  with  only  slight  modifications  to  calculate  the  wide  band  frequency  response 
with  reduced  computation  time.  A  compact  analytic  approximation  for  the  spectrum  is  obtained, 
from  which  a  direct  time  domain  (pulse)  solution  of  a  radiation  or  scattering  problem  is  easily 
derived.  Since  it  is  assumed  that  the  moment  method  model  of  the  object  has  already  been 
developed  and  tested,  little  additional  coding  work  is  required. 

As  an  example  of  the  application  of  the  hybrid  SEM/MoM  technique  presented  here,  the  radar 
cross  section  of  a  microstrlp  patch  is  calculated  over  the  1-18  GHz  band.  The  numerical 
evaluation  of  the  radar  cross  section  of  microstrlp  patches  in  the  frequency  domain  has  been 
studied  by  several  authors  14) -(5|.  The  agreement  between  measured  and  calculated  results  that 
these  authors  have  attained  has  been  very  good.  Nevertheless,  the  computationed  effort  required 
in  some  cases  has  been  quite  large.  Therefore,  several  methods  have  been  developed  to  help 
reduce  the  computational  burden.  These  methods  have  often  Involved  sparse  frequency  sampling 
or  "smart"  interpolation  schemes  based  on  additional  physical  insight  or  mathematical 
characteristics.  Interpolation  of  the  moment  method  matrix  14)  is  one  of  the  earlier  methods  that 
was  applied  with  success.  In  this  case,  the  time  required  to  fill  the  moment  method  matrix 
dominated  the  calculation  so  the  matrix  was  filled  at  only  widely  spaced  frequencies.  The 
moment  method  matrix  was  then  Interpxjlated  at  intermediate  frequencies  and  the  resulting 
matrix  equation  solved  for  the  current  at  these  intermediate  frequencies.  This  interpolation 
approach  is  not  successful  in  all  cases,  however.  The  addition  of  a  probe  attachment  basis 
function  causes  the  method  to  fail  15).  In  addition,  matrix  interpolation  schemes  do  not  help  MoM 
techniques  that  have  relatively  short  matrix  fill  times  and  longer  matrix  solution  times  (61-|7].  In 
these  situations,  a  result  of  the  matrix  equation  such  as  the  induced  current  density  rather  than 
the  matrix  equation  itself  must  be  interpolated.  However,  due  to  the  highly  resonant  nature  of 
these  quantities,  this  interpolation  is  exceedingly  dlflicult. 
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The  sections  that  follow  outline  the  basic  formulation,  discuss  the  computational 
implementation  and  present  computed  radar  cross  section  results  for  a  rectangular  air  loaded 
mlcrostrlp  patch  over  a  18:1  frequency  range. 

2.  FORMULATION  OP  THE  APPROACH 


The  details  of  the  general  approach  are  outlined  in  this  section. 

The  currents  on  a  three-dimensional  structure  such  as  the  mlcrostrlp  patch  shown  in  Figure 
1  can  be  expanded  in  terms  of  its  natural  resonances  and  natural  modes  natural 

resonances  and  corresponding  modes  are  solutions  to  the  homogeneous  problem  and  thus  are 
independent  of  any  incident  fields  or  forcing  function.  The  natural  resonances  occur  as  p>airs  in 
quadrants  1  and  11  of  the  complex  frequency  plane  and  are  s)mmetrlcally  placed  about  the 
imaginary  axis,  (when  a  frequency  plane  is  chosen  such  that  the  real  frequency  axis  corresponds 
to  continuous  undamped  waves).  The  poles  He  in  the  upper  half-plane  so  the  impulse  response 
decays  with  time. 

The  Singularity  Expansion  Method  (SEM)  current  approximation  of  the  true  spectrum  J  (r'.w) 
of  the  impulse  response  in  terms  of  the  natural  resonances  and  modes  takes  the  form 

J(r',co)  =  JsEM  = 


m«l 


1 

“-“m 


(1) 


where  are  the  mode  amplitudes  or  coupling  coefficients.  The  ampUtudes  are  determined  by 
the  excitation.  For  the  scattering  problem,  for  example,  the  amplitudes  determine  how  strongly  a 
plane  wave  incident  from  a  given  direction  excites  the  m-th  natural  mode.  TTie  natural 
resonances,  the  corresponding  natural  modes,  and  the  mode  amplitudes  for  a  particular 
excitation  are,  in  general,  unknown  and  must  be  determined  to  complete  the  current 
approximation  Jsem  (*' 

In  theory,  the  series  in  Ek^.  (1)  contains  an  infinite  number  of  natural  modes,  M  ^  and 
there  is  also  a  possible  contribution  from  an  er+lre  function.  However,  in  practice  a  small 
number  of  modes  often  suffices. 

Note  the  convenient  separation  of  the  frequency  and  incident  angle  dependencies  in  Eq.  (1). 
This  facilitates  computing  the  time  domciin  response  of  the  structure  since  once  the  currents  are 
written  in  this  form  the  multiplication  with  the  spectrum  of  the  desired  signal  and  the  inverse 
Fourier  transform  to  obtain  the  time  response  become  trivially  simple. 

The  natural  frequencies  and  modes  are  obtained  from  the  MoM  matrix  equation  that  normally 
takes  the  form 

M(r',r,a))  J{r',(o)  =  V(r,(o)  (2) 
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Figure  1.  Geometry  of  Air  Loaded  Mlcrostrlp  Patch  and  incident  Plane  Wave. 
Patch  dimensions:  a  =  36.6  mm,  b  =  26.0  mm 
Patch  height  above  ground  plane:  h  =  1.58  mm 
Plane  wave  incident  angle:  0  =  60°,  (p  =  45°. 
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where  J(r'.(o)  Is  the  current  density  vector,  M(r.r',a))  is  the  moment  method  impedance  matrix  and 
the  vector  V(r.a))  is  related  to  the  particular  forcing  function  that  may  be.  for  example,  an  incident 
plane  wave  or  voltage  gap  generator. 

The  natural  resonances  (Om  are  defined  by  the  matrix  M(r.r'.Q))  in  that  they  satisfy 

det  [M(r'.r,Q)in  )]  =  0  (3) 


which  is  a  necessauy  condition  for  a  source-free  solution.  The  correspxjnding  natural  modes 
Jhitr')  are  obtained  as  the  solutions  of  the  matrix  equation 

M(r' .  r.  (Om  ) Jm  Ir' )  =  0  (4) 

The  evaluation  of  the  unknown  amplitudes,  the  of  Ekj.  (1).  is  usually  the  major  difficulty, 
since  they  essentially  involve  evaluating  the  residue  of  M' '  V  at  to  =  cojjj.  However,  in  the  present 
case  the  available  moment  method  code  allows  a  much  simpler,  alternate  approach.  Using  this 
code  to  generate  a  CW  reference  current  (tOn)  on  the  patch  at  a  particular  real  frequency  tOn, 
the  amplitudes  can  be  determined  by  the  condition  that  the  current  estimate  JgEM  ^“n^  ^‘l- 

provides  a  least  mean  square  match  to  (t^).  that  is,  the  error 

=  PsEM  ^^^n  J  ~  "^cw  (®n  ^ 

is  minimized.  In  principle  all  the  amplitudes  determined  by  matching  at  one  single 

frequency.  However,  a  more  uniform  matching  over  the  band  of  interest  obtained  by  minimizing  a 
combined  error  from  several  frequencies 

e  =  X  ^n(Tli . )•  (6) 

n 


usually  leads  to  a  better  result. 

Once  the  currents  on  the  structure  are  known,  the  input  impedance  or  radar  cross  section 
may  easily  be  calculated. 

3.  EXAMPLE:  THE  RADAR  CROSS  SECTION  OF  A  MICROSTRIP  PATCH 

In  this  section  the  computational  issues  involved  in  applying  the  singularity  expansion 
method  to  the  study  of  scattering  from  a  microstrip  patch  will  be  discussed,  since  they  have  a 
strong  influence  on  the  overall  success  or  failure  of  the  method.  In  addition,  important  features 
of  the  method  and  its  implementation  will  be  further  outlined. 
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The  natural  resonances  were  first  determined  by  modifying  an  existing  moment  method 
computer  code  |6)-|7|  to  accept  complex  frequencies.  This  moment  method,  code  is  a  so-called 
“full-wave"  code  that  can  accurately  describe  higher-order  hybrid  current  distributions  on  the 
p>atch.  The  method  is  based  on  a  mixed  potential  integral  equation  and  expands  the  patch 
currents  in  a  series  of  "roof-top"  shaped  basis  functions.  Typically  130  basis  functions  are 
required  to  accurately  model  the  patch  currents  at  the  fourth  or  fifth  resonant  frequency.  While 
the  matrix  fill  time  is  relatively  short,  the  MoM  matrix  is  fairly  large  and  the  solution  time  of  the 
matrix  equation  grows  roughly  as  as  more  basis  functions  are  added  to  model  finer  patch 
details  or  to  study  the  antenna  at  higher  frequencies. 

The  natural  resonances  were  determined  from  Ekj.  (3),  which  implies  that  an  eigenvalue  of  the 
impedance  matrix  goes  to  zero  or.  equivalently,  the  condition  number  goes  to  infinity.  Here  the 
complex  frequency  plane  was  searched  for  maxima  of  the  condition  number  of  the  moment 
method  Impedance  matrix.  The  Muller  method  |8I  was  used  for  this  search  while  UNPACK  (9| 
routines  were  used  to  estimate  the  condition  number.  The  real  resonant  frequencies  determined 
by  the  cavity  model  were  used  as  starting  locations  for  the  Muller  method  search  since  the 
efficiency  of  this  method  depends  heavily  on  the  quality  of  the  initial  starting  locations.  The 
cavity  model  resonances  (TMpqo)  for  a  rectangular  mlcrostrlp  antenna  are  given  approximately  in 
GHz  by 


/=  0.96 


150 

ITT 


(7) 


where  a  and  b  are  the  dimensions  of  the  patch  in  millimeters.  Cj-  is  the  relative  p>ermlttlvlty  of  the 
substrate,  and  the  factor  0.96  accounts  for  the  difference  between  a  closed  cavity  and  the  actual 
microstrip .  which  has  a  fringing  field. 

The  patch  shown  in  Figure  1  was  subdivided  into  a  large  number  of  cells  (13  by  11  cells 
leading  to  a  total  of  262  rooftop  basis  functions:  132  x-directed,  130  y-directed)  in  the  moment 
method  code  to  yield  accurate  results  at  the  higher  frequencies.  The  search  for  each  complex 
resonant  frequency  required  approximately  12  evaluations  of  the  moment  method  matrix 
condition  number.  Therefore,  finding  the  first  seven  natural  resonances  to  cover  the  1-12  GHz 
band  required  approximately  84  evaluations  of  the  MoM  matrix  condition  number. 

The  natural  modes  corresponding  to  each  of  the  natural  resonances  were  found  by  use  of  a 
matrix  deflation  technique.  At  a  complex  resonance  the  natural  mode  is  a  non-trlvial  solution  to 
the  homogeneous  matrix  equation  (4).  To  determine  this  mode  one  arbitrary  mode  component  is 
set  to  a  non-zero  value,  which  leads  to  a  reduced  set  of  N-1  inhomogenous  equations  that  can  be 
solved.  This  determines  the  natural  mode  apart  from  its  magnitude,  which  can  be  set  by 
normalization.  The  numerical  condition  number  of  the  reduced  matrix  was  typically  50  times 
better  than  the  original  "numerically  singular"  matrix.  The  computation  of  each  of  the  natural 
modes,  therefore,  requires  roughly  one  evaluation  of  the  MoM  matrix  equation.  (This  method  can 
fail  if  the  reduced  matrix  equation  is  also  numerically  singular.) 
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A  second  method  for  determining  the  natural  modes  is  to  compute  the  eigenvector  associated 
with  the  minimum  eigenvalue  of  the  matrix  at  a  natural  resonance.  This  method  produces 
equally  good  results  but  is  significantly  less  efficient  computationally. 

The  same  moment  method  code  was  then  used  to  calculate  the  reference  current  distributions 
at  several  real  frequencies.  The  error  between  the  reference  moment  method  currents  and  the 
current  approximation  (Ek).  (1))  was  minimized  to  determine  the  mode  amplitudes. 

Finally,  for  the  microstrip  patch  scattering  problem  considered  here,  the  radar  cross  section 
in  the  polar  direction  {0,  «p)  was  computed  using  the  natural  mode  current  expansion  by 


®7E  =  ~  l^*^o  ®  (^x  <t>  -  Fy  COS  (}»)  /  Ejncf 

jt  ^ 

(8a) 

-  -!-|Zok^Ah  cos^0(Fx  cos(t>  +  Fy  sln(l>)/Eincf 

Jt 

(8b) 

where  Zq  Is  the  free  space  wave  Impedance,  is  the  vravenumber,  A  is  the  area  of  a  charge  cell 
on  the  patch  (that  is,  one-half  the  area  covered  by  a  roof-top  basis  function),  h  is  the  height  above 
the  ground  plane,  and  and  Fy  are  array  factors  for  the  rooftop  basis  functions  on  the  patch 
given  by 


f'x  =  X  (9a) 

m 

=  X  (9b) 

n 

The  currents  and  Jy^  are  the  x  and  y  directed  components  of  the  current  approximation 
(1),  and  V  is  the  angle  between  the  vectors  to  the  field  point  r  and  to  the  source  p>oint  r'.  In  this 
derivation  it  is  assumed  that  h  «  X  so  that  the  array  factor  for  the  patch,  and  its  negative  image 
in  the  groundplane,  is 

e-jk^hcose  _  eJk^hcosG  -  .j2kohcose 

The  far  field  for  a  patch  with  known  current  is  then  calculated  following,  for  instance,  van 
Bladel  [10]. 

The  computational  effort  to  evaluate  the  mode  amplitudes  was  negligible  while  the 
computation  of  the  radar  cross  section  of  the  patch  using  the  current  estimates  over  the  entire 
band  was,  for  the  case  chosen  here,  roughly  equivalent  to  one  evaluation  of  the  MoM  matrix 
equation.  Thus,  with  the  present  state  of  affairs,  the  computational  effort  required  to  Implement 
the  SEM  and  MoM  solutions  was  roughly  equivalent  at  approximately  90  frequencies.  SEM  is 
more  efficient  for  a  larger  number  of  frequency  points  while  MoM  is  better  for  a  smaller  number  of 
frequencies. 
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Both  the  moment  method  and  the  singularity  expansion  method  can  compute  radar  cross 
section  efOclently  for  multiple  Incident  angles.  The  moment  method  approach  must  Invert  the 
MoM  matrix  once,  compute  the  multiple  rlght-hand-sldes  and  multiply.  The  present  technique 
requires  the  MoM  solution  for  the  reference  current  at  several  frequencies  for  each  incident  angle 
to  determine  the  mode  amplitudes. 

4.  RESULTS  AND  DISCUSSION 

The  radar  cross  section  of  the  air-loaded  patch  of  Figure  1,  computed  by  using  the  natural 
mode  current  expansion.  Is  presented  in  this  section  and  compared  to  the  direct  calculation  by 
the  moment  method.  While  In  principle  the  present  method  and  the  MoM  code  can  be  applied  to 
dielectric  loaded  patch  antennas  where  the  substrate  of  the  antenna  is  not  air,  questions  about 
surface  wave  pole  migration  in  the  complex  plane  when  studying  the  antenna  at  complex 
frequencies  were  beyond  the  scope  of  the  current  work.  The  precise  location  of  the  surface  wave 
poles  must  be  known  during  the  calculation  of  the  Sommerfeld  integrals  that  are  the  Green's 
functions  of  the  integral  equation. 

Figure  2  shows  the  locations  of  the  first  11  members  of  a  family  of  complex  natural 
resonances  that  lie  near  the  real  frequency  axis.  The  real  resonant  frequencies  predicted  by  the 
cavity  model  that  were  used  as  starting  locations  for  the  Muller  method  search  and  the 
corresponding  mode  numbers  are  also  shown.  Note  the  surprisingly  close  correspwndence 
between  the  real  and  complex  resonances.  (The  correspondence  was  established  by  comparing  the 
mode  current  distributions).  A  list  of  15  poles  that  were  found  In  this  family  is  given  in  Table  1. 

A  second  family  of  natural  resonances,  shown  in  Figure  3,  lies  in  the  complex  plane  far  above 
the  first.  These  modes  play  a  secondary  role  in  the  late  time  impulse  response  of  the  patch  and 
were  more  difficult  to  find  since  there  exists  no  clear  starting  point  for  the  search.  The  second 
family  modes  were  not  Included  in  the  results  that  follow. 

The  convergence  of  the  first  family  of  poles  with  the  number  of  moment  method  basis 
functions  is  shown  in  Figure  4.  The  subdivision  of  the  patch  into  a  larger  number  of  cells 
improves  the  convergence  of  the  natural  modes  at  the  higher  frequencies  but  does  not 
significantly  improve  the  final  radar  cross  section  results.  The  number  of  cells  chosen  here,  13  x 
11,  corresponds  to  about  9  cells  per  wavelength  at  12  GHz  and  6  per  wavelength  at  18  GHz. 

Figure  5a  shows  a  quiver  diagram  of  the  in-phase  and  quadrature-phase  components  of  the 
seventh  natural  mode.  The  complex  frequency  of  this  mode  is  approximately  1 1.20  +  J0.33  GHz 
and  corresponds  to  the  TM220  cavity  mode  resonant  at  f  =  1 1.52  GHz.  We  note  that  the  in-phase 
component  is  roughly  8  times  stronger  than  the  quadrature  component.  Figure  5b  similarly 
shows  the  15:th  natural  mode,  which  clearly  is  seen  to  corresponds  to  the  TM500  cavity  mode. 
This  time  the  in-phase  component  dominates  by  a  factor  of  1 1. 
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Real  Frequency  (GHz) 

•  Natural  Resonartcas 

o  Cavity  Resonances 


Figure  2.  First  Family  Pole  Locations  in  the  Complex  Frequency  Plane  for  the  Patch  of  Figure  1. 

The  corresponding  resonant  frequencies  predicted  by  the  cavity  model  are  shown  along 
the  real  axis. 
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Figure  3.  Second  Family  Pole  Locations  in  the  Complex  Frequency  Plane  for  the  Patch  of  Figure  1 
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Figure  4.  Convergence  of  the  First  Family  of  Poles  In  the  Complex  Frequency  Plane  for  Different 
Numbers  of  Moment  Method  Basis  Functions. 
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Figure  5a.  In-phase  and  Quadrature  Phase  Current  Components  for  the  Seventh  Natural  Mode 
of  the  Patch  Shown  in  Figure  1.  Corresponding  to  the  TM120  Cavity  Mode. 


12 


Max.  Value  =  0.799 x10**-1  Amps 


Figure  5b.  In-phase  and  Quadrature  Phase  Current  Components  for  the  15;th  Natural  Mode 
Corresponding  to  the  TM500  Cavity  Mode. 
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Figures  6  through  7  compare  the  SEM  and  MoM  radar  cross  section  (000)  predictions  for  the 

air  loaded  patch  of  Figure  1  using  seven  pmles  and  natural  modes.  The  RCS  was  calculated  at 
101  frequencies  in  all  cases  to  produce  a  fairly  smooth  curve.  Note  that  as  more  matching  points 
are  added  the  quality  of  the  radar  cross  section  prediction  improves.  This  implies  that  a  slight 
decrease  in  approximation  accuracy  at  one  or  two  frequencies  is  traded  for  a  greatly  improved 
average  accuracy  over  a  wider  band. 


Table  1.  Cavity  Resonances  and  Corresponding  Complex  Poles 


Pole 

number 

Mode 
indices 
(p.  q.  0) 

Cavity 

resonance 

(GHz) 

Complex  Pole 
(GHz) 

1 

100 

3.93 

3.811+j0.515 

2 

010 

5.54 

5.207+J0. 1463 

3 

110 

6.79 

6.532+j0.0997 

4 

200 

7.87 

7.63  IfjO.  1658 

5 

210 

9.62 

9.343+j0.1488 

6 

020 

11.1 

10.506+J0.3794 

7 

120 

11.75 

11.200+j0.3324 

- 

300 

11.80 

? 

8 

310 

13.0 

12.636+j0.2323 

9 

220 

13.6 

13.087+j0.3125 

10 

400 

15.7 

15.079+j0.3551 

11 

320 

16.2 

15.634+j0.3403 

12 

410 

16.7 

16.049+j0.3141 

13 

030 

16.6 

16.2 19+jO. 5028 

- 

130 

17.1 

? 

14 

230 

18.4 

17.566+J0.4847 

- 

420 

19.2 

? 

15 

500 

19.7 

18.561+J0.3965 

A  detailed  picture  of  the  quality  of  the  current  approximation  obtainable  with  seven  modes  is 
given  in  Figure  8  which  compares  the  patch  current  Jsem  'true'  patch  current 

Jmom  computed  by  the  MoM.  The  relative  error  was 


£  =  l*%EM  ^l'^o.24. 
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RCS(dBsc|ni) 


Figure  6.  Comparison  of  the  RCS  Calculated  Using  the  SEM  Current  Model  and  the  Exact  MoM 
Model  for  the  Patch  of  Figure  1.  The  SEM  current  model  Is  matched  to  the  MoM 
currents  at  one  frequency  (4  GHz  /  curve  SEMI)  and  two  frequencies  (4,  6  GHz  /  curve 
SEM2).  Incident  and  scattered  fields  TM  to  the  z-axls. 
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RCS(dBs<|n} 


Figure  7.  Comparison  of  the  RCS  Calculated  Using  the  SEM  Current  Model  and  the  Exact  MoM 
Model  for  the  Patch  of  Figure  1.  The  SEM  current  model  is  matched  to  the  MoM 
currents  at  four  frequencies  (4.  6,  8,  10  GHz  /  curve  SEM4).  Incident  and  scattered 
fields  TM  to  the  z-axis. 
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SEM 


Figure  8.  In-phase  Components  of  Patch  Current  as  Obtained  With  Exact  MoM  ComputaOon  fTop) 
and  With  a  Seven  Mode  Expansion  (Bottom),  at  10  GHz.  Incident  field  TM  to  the  z-axis. 
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frequency  (GHi) 


Figure  9.  Comparison  of  the  RCS  Calculated  Using  the  SEM  Current  Model  and  the  Exact  MoM 

Model  for  the  Patch  of  Figure  1.  The  SEM  current  model  includes  15  natural  modes  and 
is  matched  to  the  MoM  currents  at  eight  frequencies  (2,  4,  6,  8,  10,  12.  14,  16  GHz). 
Incident  and  scattered  fields  TM  (top)  and  TE  (bottom)  with  respect  to  the  z-axls. 


Although  in  this  case  we  matched  at  one  frequency  only,  this  is  also  a  typical  value  when 
matching  at  several  frequencies. 

It  is  noteworthy  that  the  high  frequency  limit  of  the  current  spectrum  can  be  readily  extended 
with  more  poles  and  natural  modes.  The  use  of  15  poles  and  modes  in  the  current  expansion 
JsEM  gives  excellent  agreement  between  the  corresponding  radar  cross  section  and  the  MoM 
computation  up  to  about  18  GHz  as  shown  in  Figure  9.  The  small  difference  observed  in  Figure  9 
(top)  between  11-12  GHz  may  be  an  effect  of  the  missing  counterpart  of  the  (3(X))-cavlty 
resonance,  which  we  were  unable  to  find  in  the  complex  plane.  Likewise  the  increasing  difference 
between  the  two  curves  at  the  high  frequency  end  is  due  to  an  Incomplete  number  of  poles  in  that 
vicinity,  (in  this  example  the  MoM  values  may  have  somewhat  reduced  accuracy  since  at  18  GHz 
only  6  basis  functions  per  wavelength  are  used.  However,  this  is  less  important  since  the  present 
purpose  is  to  show  only  that  the  natural  mode  expansion  does  approximate  the  MoM  curve  wellj. 

In  the  low  frequency  regime  below  the  first  resonance,  the  radar  cross  section  of  the  patch 
varies  theoretically  as  k®  (24  dB  per  octave),  as  shown  in  Appendix  A.  The  MoM  based  curve 
follows  this  behavior.  The  natural  mode  expansion  as  described  by  Eq.  (1)  however,  does  not, 
naturally  have  the  correct  low  frequency  behavior,  in  fact  it  does  not  even  go  to  zero  when  the 
frequency  is  zero.  By  expanding  the  array  factor  of  Eq.  (9)  in  a  Taylor  series  about  k^  =  0  it  can  be 
shown  that  the  correct  ko  -dependence  requires 

3  «%EM («•'.<»>  =  0)_Q  jjQj 

do) 

but,  unfortunately,  these  additional  conditions  are  not  readily  imposed  on  Jsem-  ^  partial 
correction  is  obtained  by  subtracting  the  DC-value  natural  mode  expansion 

(1).  This  leads  to  a  k®  -dependence  but  does  not  seem  to  significantly  improve  the  overall 
approximation  accuracy,  as  shown  in  Figure  10. 

Finally,  we  note  that  once  the  patch  current  is  available  the  radiation  pattern  or  blstatlc  radar 
cross  section  for  any  desired  direction  is  easily  computed.  An  example  of  the  blstatlc  RCS  of  the 
patch  at  16  GHz  is  given  in  Figure  11.  This  frequency  is  sufficiently  high  for  the  scattering  pattern 
to  exhibit  several  lobes  and  these  are  apparently  quite  well  reproduced  by  the  natural  mode 
expansion. 

5.  CONCLUSION 

A  method  that  uses  existing  frequency  domain  MoM  codes  to  obtain  a  wideband 
characterization  of  three-dimensional  antennas  or  scatterers  has  been  developed.  The  codes  must 
be  extended  to  complex  frequencies,  although,  in  general,  this  is  a  simple  matter.  The  resulting 
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RCS(dBs^) 


Figure  10.  Same  Case  as  Figure  7,  Except  That  the  SEM  Current  Model  Has  Been  Modifled  to 
Include  a  DC  Correction  Term. 
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BISTATIC  RCS  (dBsqn) 


8ISTATIC  RCS  (dBscyn) 


r  igure  11.  BIstatic  RCS  at  16  GHz  Calculated  Using  the  SEM  Current  Model  and  the  Exact  MoM 
Model  for  the  Patch  in  Figure  1.  The  scattering  angle  is  measured  in  the  plane  of 
Incidence,  with  30°  being  the  direction  of  specular  relectlon.  TM-TM  and  TM-TE  denote 
CO-  and  cross-polarized  RCS.  resp. 
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current  model  Is  very  concise  and  numerically  efficient.  For  example,  we  used  only  7  current 
vectors  to  represent  the  complex  spectrum  shown  In  Figure  7.  as  compared  to  the  100  vectors 
required  by  conventional  dense  sampling. 

The  natural  frequencies  and  modes  provide  a  complete  analytical  model  from  which  the 
complex  time  domain  responses  are  readily  obtained  over  all  angles.  A  change  In  incidence  angle 
requires  a  recomputation  of  the  mode  amplitudes  only,  which  Is  a  minimal  computational  effort. 

In  the  example  considered  we  obtained  an  accurate  analytic  representation  for  the  spectrum 
of  the  patch  current  up  to  18  GHz.  which  Is  commensurate  with  a  spatial  pulse  width 
Ax  =  c  At  —  c  /  fjnax  ~  about  half  the  scatterer  length  (3.66  cm).  Thus  we  are 

begliming  to  actually  resolve  the  scatterer  spatially  and  are  doing  considerably  better  than  a  'late 
time’  response,  valid  only  after  the  pulse  has  painted  over  the  entire  object. 

The  method  has  been  shown  to  work  well  for  three  dimensional  structures,  which  are  highly 
resonant,  l.e..  with  poles  close  to  the  real  axis.  It  remains  to  be  seen  what  the  effects  of  a  lower  Q 
would  be. 

At  this  point  it  is  premature  to  compare  the  computer  times  of  the  proposed  approach,  based 
on  a  natural  mode  expansion,  with  a  straightforward  MoM  approach,  using  dense  frequency 
sampling.  In  our  approach  most  of  the  computer  time  was  expended  in  searching  the  complex 
plane  for  the  natural  resonances  for  which  we  used  standard  routines  from  a  computer  llbraiy.  In 
contrast,  the  MoM  codes  have  been  highly  optimized  for  computational  efficiency.  Thus  the  key  to 
the  overall  success  of  the  proposed  method  is  an  efficient  root  finding  routine  and  there  Is  both 
the  need  and  the  room  for  significant  Improvement  in  this  respect. 
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Appendix 

The  ko-dependence  of  the  Radar  Cross  Section  of  a 
Microstrip  Patch  at  Low  Frequencies 

The  scattering  from  a  metallic  plate  over  a  reflecting  ground  plane  can  be  anal3rzed  in  terms  of 
scattering  from  two  plates  in  free  space  that  are  Illuminated  from  above  and  below.  In  the  upper 
half-space  the  total  fields  are  identical  for  the  original  and  equivalent  problems  shown  in  Figure 
Al. 

First  consider  a  single  plane  wave  incident  on  the  equivalent  structure  from  above.  For  low 
frequencies  the  current  on  the  upper  patch  can  be  expanded  as 

J*  =  Jo  +koJ;  J2  +.  .  .  (Ala) 

and  on  the  lower  patch 

j"=j; +koj; +k:2  j;  + . . .  (Aib) 

where  kg  is  the  wavenumber.  The  first  terms  in  each  expansion  are  the  same  since  both 
patches  have  the  same  currents  as  kg  0. 

For  incidence  from  below,  due  to  s)anmetiy  and  the  sign  reversal  of  the  incident  field,  the 
upper  patch  will  carry  the  current  -  J"  and  the  lower  -  J*.  Superimposing  these  two  cases  leads 
to  upper  and  lower  patch  currents  AJ  and  -AJ,  respectively,  where 

AJ  =  ko(j;  -j;  ).  (A2) 


The  vector  potential  due  to  currents  on  the  equivalent  structure  is 


A 


=  iiQf 

4«  JSj  +  S2 


J 


(A3) 


where  the  integration  is  over  both  patches  Sj  and  S2.  Substituting  the  p>atch  currents  above. 
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using  the  fact  that  the  patch  dimensions  are  small  relative  to  the  wavelength,  and  making  the 
usual  far  field  af^nudmatlons  leads  to 


A  =  -Hk  (-J2kohcose  )  L  AJ  dS.  (A4) 

4n  r  JSi 

We  must  now  determine  the  kg-dependence  of  the  integral  appearing  in  Ek).  (A4). 

For  a  single  incident  plane  wave  the  patch  currents  are  equal  to  first  order  and  thus  the  far- 
fleld  vector  potential  at  low  frequencies  becomes 


Jo  dS  (A5) 

4jc  r  JSi  u 

The  electric  far  field  is  proportional  to  kg  times  the  vector  potential 

E<p  “  ^ 

Ee  “  ^0^0  (A6) 

and  from  Rayleigh  scattering  theory  we  know 

oc  kl  (A7) 

This  leads  to  the  result  that  the  integral  appearing  in  Eq.  (A5)  is  proportional  to 

Jg^JodSocko  (A8) 

Since  the  integrals  over  Jq  and  ,  Jj  are  of  the  same  order  as  kp  the  integral  in  Eq.  (A4)  is 
proportional  to  kp  . 

Jg^  AJ  dS  =  kp  Jg^  -j{)dS^kl  (A9) 

Substituting  Eq.  (A9)  in  Eq.  (A4)  and  using  Eq.  (AS)  leads  to  a  far  field  for  the  patch  over  a 
ground  plane 
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-<p 


*^0 


(10) 


and  a  low  frequency  radar  cross  section  proportional  to  Ic* 

o(0.<p) «  Ic®  . 


(All) 


This  corresponds  to  a  roll-off  of  24  dB  per  octave. 


Original  Problam 


Equivalent  Problem 


Figure  Al.  Scattering  From  a  Microstrip  Patch  Related  to  a  Patch  Pair  in  Free  Space. 
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Rom«  Laboratory  plans  and  executes  an  interdisciplinary 
progrjMi  in  research,  development,  test,  and  technology 
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Commtmications  and  Intelligence  (C3I)  activities  for  all 
Air  Force  platforms.  It  also  executes  selected 
acquisition  programs  in  several  areas  of  expertise. 
Technical  and  engineering  support  within  areas  of 
competence  is  provided  to  ESC  Program  Offices  (POs)  and 
other  ESC  elements  to  perform  effective  acquisition  of 
C3I  systems.  In  addition,  Rome  Laboratory's  technology 
supports  other  AFMC  Product  Divisions,  the  Air  Force  user 
community,  and  other  OOD  and  non>DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research 
programs  in  areas  including,  but  not  limited  to, 
communications,  command  and  control,  battle  management, 
intelligence  information  processing,  computational 
sciences  and  software  producibility,  wide  area 
surveillance/ sensors,  signal  processing,  solid  state 
sciences,  photonics,  electromagnetic  technology, 
superconductivity,  and  electronic 
reliability/maintainability  and  testability. 


